Solving second order ordinary differential equations by 

extending the PS method 

L.G.S. Duartet L.A. da Mota^and J.E.F.Skea* 
February 4, 2008 

Abstract 

An extension of the ideas of the Prehe-Singer procedure to second order differ- 
ential equations is proposed. As in the original PS procedure, this version of our 
method deals with differential equations of the form y" = M{x,y,y')/N{x,y,y'), 
where M and N are polynomials with coefficients in the field of complex numbers C. 
The key to our approach is to focus not on the final solution but on the first-order 
invariants of the equation. Our method is an attempt to address algorithmically the 
solution of SOODEs whose first integrals are elementary functions of x, y and y'. 
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1 Introduction 



The fundamental position of differential equations (DEs) in scientific progress has, 
over the last three centuries, led to a vigorous search for methods to solve them. 
The overwhelming majority of these methods are based on classification of the DE 
into types for which a method of solution is known, which has resulted in a gamut 
of methods that deal with specific classes of DEs. This scene changed somewhat 
at the end of the 19th century when Sophus Lie developed a general method to 
solve (or at least reduce the order of) ordinary differential equations (ODEs) given 
their symmetry transformations [Q, ||, |3|. Lie's method is very powerful and highly 
general, but first requires that we find the symmetries of the differential equation, 
which may not be easy to do. Search methods have been developed Q ^ to extract 
the symmetries of a given ODE, however these methods are heuristic and cannot 
guarantee that, if symmetries exist, they will be found. 

On the other hand in 1983 Prelle and Singer (PS) presented a deductive method 
for solving first order ODEs (FOODE) that presents a solution in terms of elemen- 
tary functions if such a solution exists Q. The attractiveness of the PS method 
lies not only in its basis on a totally different theoretical point of view but, also in 
the fact that, if the given FOODE has a solution in terms of elementary functions, 
the method guarantees that this solution will be found (though, in principle it can 
admittedly take an infinite amount of time to do so). The original PS method 
was built around a system of two autonomous FOODEs of the form x = P{x,y), 
y = Q{x,y) with P and Q in C[x,y] or, equivalently, the form y' = R{x,y), with 
R{x, y) sl rational function of its arguments. Here we propose a generalization that 
allows us to apply the techniques developed by Prelle and Singer to second order 
differential equations (SOODEs). The key idea is to focus not on the final solution 
of the equation, but rather its invariants. 

This paper is organized as follows: in section the reader is introduced to the PS 
procedure; section ^ addresses our approach extending the ideas of the PS procedure 
to the case of SOODEs and discusses how generally applicable the method is to such 
equations. Section ^ is dedicated to some examples solved via our procedure and, 
finally, conclusions are presented in section 0. 



2 The Prelle-Singer Procedure 

Despite its usefulness in solving FOODEs, the Prelle-Singer procedure is not very 
well known outside mathematical circles, and so we present a brief overview of the 
main ideas of the procedure. 

Consider the class of FOODEs which can be written as 

, ^dy ^ M{x,y) 
^ dx N{x,y) ^ ' 

where M(x, y) and N{x, y) are polynomials with coefficients in the complex field C . 

In 1^], Prelle and Singer proved that, if an elementary first integral of (0) exists, 
it is possible to find an integrating factor R with i?" G C for some (possible 
non-integer) n, such that 

^ + ^=0. (2) 

ax oy 
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The ODE can then be solved by quadrature. From (^) we see that 

^BR ^dN ^BR ^dM ^ 
ox ox oy oy 



+ (4) 



Thus 

D[R] _ fdN dM' 

R \ dx dy 
where 

Now let R = Y\i where fi are irreducible polynomials and rij are non-zero 
integers. From (^), we have 

D[R] _ _E^/r"'™^^[/dn,v./r' 



R n./r Ukfr 

From (Q), plus the fact that M and N are polynomials, we conclude that D\R\IR 
is a polynomial. Therefore, from (|6|), we see that 

We now have a criterion for choosing the possible fi (build all the possible 
divisors of and, by using (Q) and (^), we have 

^ fi \dx^ dy )■ ^ ' 

If we manage to solve (0) and thereby find n^, we know the integrating factor for 
the FOODE and the problem is reduced to a quadrature. Risch's algorithm can 
then be applied to this quadrature to determine whether a solution exists in terms 
of elementary functions. 



3 Extending the Prelle-Singer Procedure 

In the previous section, the main ideas and concepts used in the Prelle-Singer pro- 
cedure were introduced. Here we present an extension of these ideas applicable to 
SOODEs. The main idea is to focus on the first order invariants of the ODE rather 
than on the solutions. 



3.1 Introduction 

Consider the SOODE 

" = ^ = Mj^^y^y') /ox 

^ dx^ N{x,y,y')' 

where M(x, y, y') and N{x, y, y') are polynomials with coefficients in C . We assume 
that (^) has a solution in terms of elementary functions, in which case there are two 
independent elementary functions of x y and y' which are constant on all solutions 
of (^), namely the first order invariants 

Ii{x,y,y') = Ci i = l,2. (9) 
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Without loss of generalization we consider one of these and, dropping the index on 
li we have 

dl dl dl 
dl = — dx + — dy + — dy' = 0. (10) 
ox oy ay' 

Now, introducing the notation ^ = /„, we have 

h + Iyy' + Iy>y" = (11) 

and so 

/ = -h±h!l^ (12) 

lyl 

which is (^ in terms of the differential invariant /. Rewriting (^) as 

M 

— dx-dy' = (13) 

and observing that 

y' dx = dy, (14) 



we can add the identically null term S{x, y, y')y' dx—S{x, y, y') dy to (|13|) and obtain 
the 1-form 

^ + Sy''] dx-Sdy-dy' = 0. (15) 



Notice that the 1-form (|16| ) must be proportional to the 1-form (|10[). So, since 
the 1-form (10) is exact, we can multiply (|l6|) by the integrating factor R{x,y,y') 
to obtain 

dl = R{(^ + Sy')dx- RSdy- Rdy' = (16) 

where cj) = M/N. 

Comparing equations (|lO|) and ([T^), 

4 = R{(l) + Sy'), 

ly = RS, 

ly, = -R. (17) 

Now equations (|T7| ) must satisfy the compatibility conditions Ixy = Iyx,Ixy' = l^y'x 
and lyy' = lyiy. Thls implies that 

D[S] = -c^y + S<l)y, +3"^, (18) 
D[R] = -R{S + <Py'), (19) 

Ry = Ryl S -\- Sy'R, (^O) 

where the differential operator D is defined as 

d f d d 

dx ~^ ^ dy ^ dy' ^ ^ 

Combining ( p^ ) and (|^) we obtain 

D[RS] = -R(^y. (22) 

So if the product of S and the integrating factor is a rational function of x, y 
and y', then D[RS] is too. Since (j) is rational (and so, therefore, \s (j)y), equation (22) 
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tells us that R is rational. Using (19) and similar arguments we conclude that S 



must be a rational function of x, y and y' . 

In summary, from ( p^ ) it follows that the supposition that RS is rational can 
be equated to the existence of a first order invariant whose derivatives in relation 
to X, y and y' are rational functions. With this in mind we restate the original 
supposition in the form of a conjecture. 

3.2 The Conjecture 

We first state a result proved in |Q: Theorem: Let K be a differential field of 

functions in n + 1 variables and L an elementary extension of K . Let f be in K 
and assume there exists a nonconstant g in L such that g is constant on all solutions 
of y^^^ = f{x,y,y',y",...,y^'^^^^). Then there exist wi, ... ,Wm algebraic over K 
and constants Ci, . . . ,Cm such that 

woix, y, y', y", .., y^""!)) + ^ c, logiw^ix, y, y' , y" , (23) 

i 

is a constant on all solutions of y*^"^ = f{x, y, y', y" , . . . , y^"~^^). 

This result shows that for the particular case of SOODEs whose solutions are 
elementary, there are two independent first order invariants of the form 

wo{x,y,y') +^CilogWi{x,y,y'). (24) 

i 

Our conjecture is that if these two first order invariants exist it is always possible 
to find a function of them (which will, therefore, itself be a first order invariant) of 
the form 

zo{x,y,y') + ^Cilog[zi{x,y,y')], (25) 

i 

where zt are rational functions x, y and y'. 

ConjecturerZei K be a differential field of functions in three variables and L an ele- 
mentary extension of K . Let f be in K and assume there exist two independent non- 
constant {^1,52} in L such that gi are constant on all solutions of y" = f{x,y,y'). 
Then there exists at least one constant of the form 

zo{x,y,y') + ^Cilog{zi{x,y,y')) (26) 

i 

where the Zi are in K . 

By the previous reasoning it can be seen that ( p6[ ) implies that the product RS 
is a rational function of x, y and y' . 

If this conjecture holds, then our extension of the PS method applies to all 
SOODEs of the form (^). Though we have not been able to prove our conjecture, 
extensive trials while developing this procedure has not revealed any counter exam- 
ple. Even if the conjecture is false, our experience with real test cases has shown 
that the method is, at least, applicable to the vast majority of SOODEs of the 
form (II. 
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3.3 Finding R and S 

Our conjecture implies that, if the SOODE to be solved has an elementary general 
solution, then S" is a rational function which we may write as 

We can also see that (|^) does not involve R. So, given a degree bound on the 
polynomials Sn and Sd, we may find a set of solutions to this equation which are 
then candidates to solve the system of equations (p^)~(|20|). 

From ([l^) we have 



D[R] (o^^ ^ (M\ SnN^ + Sd{NMy.-MNy,) 



which can be rewritten as 
V[R] 



-SnN^ + Sd{NMy, - MNy,), (29) 



R 

where the differential operator D is defined as 

V = {SdN^)D. (30) 

We keep in mind that 

• Sn, Sd, N and M are polynomials in x, y and y'; 

dx ' dy dy 



P is a linear differential operator whose coefficients of j?- and -S-r are 



polynomials in x, y and , y'; 
• i? is a rational function of x, y and y', which we may write as 



Rn ^ T,i,j,kCijkx'y^y'- 
T,i,j,k dijkx^y^y'' 



R = t^ = ^M2!iZ!d^_JLJL__ (31) 



If we have a theoretical limit on the degrees of Rm and Rd (a degree bound), we 
may use a procedure analogous to that described in section |2| to obtain candidates 
for the integrating factor R. We simply construct all polynomials in x, y and y' up 
to the degree bound. 



3.4 Reduction of the SOODE 

Once R and S have been determined using equations ( p!7[ ) we have all the partial 
first derivatives of the first order differential invariant, I{x,y,y'), which is constant 
on the solutions. This invariant can then be obtained as 

I{x,y,y') = JR{(P + Sy')dx - 
J[RS + ^Jr{(P + Sy') dx]dy - 
^(^lR{<P + Sy')dx- J[RS+-^ lR{(t> + Sy')dx]dy^ dy' . (32) 



R + 
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The equation I(x,y,y') = Ci can then be solved to obtain a FOODE for y': the 
reduced ODE 

y' = ip{x,y,Cl). (33) 

To obtain the general solution of the original ODE, we can apply the Prelle-Singer 
method in its original form to this reduced ODE. Thus, if our conjecture is correct, 
the method proposed here (for SOODEs of the form (^)) is as algorithmic as the 
original PS method for FOODEs. We note that the original PS method fails to 
be what is strictly an algorithm because no theoretical degree bound is yet known 
for the candidate polynomials which enter in the prospective solution, and so the 
procedure has no effective terminating condition for the case when an elementary 
solution does not exists. In practice, a terminating condition is put in by hand 
(it is found that polynomials of degree higher than 4 lead to computations which 
are overly complex for the average desktop computer). However, should such a 
degree bound be established, and our conjecture shown to be true, then the method 
proposed here would be an algorithm for deciding whether elementary solutions of 
SOODEs of the form (g) exist. 



4 Examples 



In this section we present examples of physically motivated SOODEs that are solved 
by our procedure]]. As a simple illustrative example, we begin with the classical 
harmonic oscillator and then consider some nonlinear SOODEs which arise from 
astrophysics and general relativity. 

Example 1: The Simple Harmonic Oscillator 
In its simplest form, the equation for the simple harmonic oscillator is 

y" = -y- (34) 

For this ODE equations (|l|), (^) and (|0|) are 

S:, + y'Sy-ySy> = l + S\ (35) 
R^ + y'Ry-yRy' = -RS, (36) 

Ry- Ry^S - Sy^R = 0. (37) 

One possible solution to these equations is 

S=^, R = y'. (38) 

y' 

From this, and using (^), we get the reduced ODE 

Cl = 2/2 + y'^ (39) 

which, of course, represents the energy conservation for the oscillator. 

This example is very simple and leads to a form of (p which is independent 
of X and y' . And, as with all linear ODEs, alternative and more straightforward 
solution methods exist. The other examples illustrate the solution method at work 
for non-linear SOODEs which can be placed in the form (S). 



^We present only the reduction of the SOODEs since the integration of the resulting FOODE can be 
achieved by various methods, including the PS method itself. 



7 



Example 2: An Exact Solution in General Relativity 

A rich source of non-linear DEs in physics are the highly non-linear equations of 
General Relativity. In general, Einstein's equations are, of course, partial DEs, but 
there exist classes of equations where the symmetry imposed reduces these equations 
to ODEs in one independent variable. One such class is that of static, spherically 
symmetric solutions for stellar models, which depend only on the radial variable, r. 
The metric for a general statically spherically spacetime has two free functions, A(r) 
and //(r) say. On imposing the condition that the fluid is a perfect fluid, Einstein's 
equations reduce to two coupled ODEs for A(r) and /x(r). Specifying one of these 
functions reduces the problem to solving an ODE (of first or second order) for the 
other. 

Following this procedure, Buchdahl |^] obtained an exact solution for a relativis- 
tic fluid sphere by considering the so-called isotropic metric 

S^ = {1- /)2(1 + f)-H^ _ (1 + /)4[^2 ^ ^202 ^ ^.^2 Q ^2] 

with / = /(r). The field equations for /(r) reduce to 

//" - 3/'2 - r-V/' = 0. 
Changing notation with y{x) = f{r), equations ([l^), (|l9|) and (|20|) assume the form 



c , /c ,y'i^y'x + y)Q y' , y' i^y' x + y) 

Jx+yjy-\ — ^y' = ^;72 



xy xy xy^ 



^y'- + y +3y-)s + s\ (40) 



xy 

Rx + y'Ry+ '' ^'''^^'^ Ry, = -Rfs+'-y^ + S^-), (41) 
xy V xy yj 

Ry- RylS - SylR = 0. (42) 

One solution of those equations is 

5 = -3-, R=^ (43) 
y xy-^ 



By using (32) we obtain the reduced FOODE: 

Ci = y'/{y^x) (44) 
which is separable and easily integrated to obtain the general solution 

y{xf = (-Cix^ + C2y\ (45) 

Example 3: A Static Gaseous General-Relativistic Fluid Sphere 

In a later paper Q, Buchdahl approaches the problem of the general relativistic 
fluid sphere using a different coordinate system from the previous example. For ease 
in comparison of the originals. 

Substituting the ^(r) of the original Writing y{x) instead of the ^(r) in the 
original, arrives at the equation 

x'^y 
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For this SOODE, eqs (0, and ^ become: 



bx+yby^ O Oyl = -2X H Ty-^ 

+2 ^ + 52 (47) 

y 

Rx + y'Ry + ''^' \^ ~^ Ry' = -r(s + 2^) (48) 

\ yj 

Ry- Ry'S - Sy'R = (49) 
One solution to those equations is: 

-x^y'^ - xyy' + 1 y + xy' 

^= 2,2/ ' ^= 2~- (^0) 

xy^ + x'^yy' xy'' 

From this, using eq. (|32|), we get the reduced FOODE: 

_2xyy' + y'^ + x'^y'^ -1 

which can be solved to: 



yix)' = tMV2^iC2+x)r 

'2 Ci +2t&n{V2^/CJ{C2+x)yCl) x"^ 



This example has an extra feature: It is not solved by other solvers we have 
tried (mainly the Maple solver, in the version 5, that we believe to be the best). So, 
apart from the (already) very interesting fact that our approach is an algorithmic 
attempt to solve SOODEs, we have also this present fact, i.e., some SOODEs are 
solved via our method and "escape" from other very powerful solvers. 



5 Conclusion 



In this paper, we presented an approach that is an extension of the ideas devel- 
oped by Prelle-Singer to tackle FOODEs. We believe it to be the first technique 
to address algorithmically the solution of SOODEs with elementary first integrals. 

Here, we dealt with a restrict class of SOODEs (namely, the ones of the form 
(P)). However, we can use our method in solving SOODEs where (f){x, y, y') depends 
on elementary functions of x,y,y', following the developments for the Prelle-Singer 



approach for FOODEs ||10|, 11|. We are presently working on those ideas. 



The generality of our approach is based on a conjecture (see section (|3.2| )) that 
we have already proved for many special cases. Even if the conjecture is proven false, 
our approach is a powerful tool in dealing with SOODEs since we have extensively 
tested it with many equations, both from mathematics and physical origin. In fact, 
since all the examples we have encounter have been solved by our approach, we are 
preparing a computational package implementing the Prelle-Singer procedure (and 
our present extension) to be submitted to Computer Physics Communications. 
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